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Abstract 

We consider a class of stochastic processes containing the classical and well-studied 
class of Squared Bessel processes. Our model, however, allows the dimension be a 
function of the time. We first give some classical results in a larger context where a 
time-varying drift term can be added. Then in the non-drifted case we extend many 
results already proven in the case of classical Bessel processes to our context. Our 
deepest result is a decomposition of the Bridge process associated to this generalized 
squared Bessel process, much similar to the much celebrated result of J. Pitman and 
M. Yor. On a more practical point of view, we give a methodology to compute the 
Laplace transform of additive functionals of our process and the associated bridge. 
This permits in particular to get directly access to the joint distribution of the value 
at t of the process and its integral. We finally give some financial applications of our 
results. 



Keywords Squared Bessel process; Bessel bridges decomposition; Laplace transform; 
Levy Ito representation; Financial applications. 
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Introduction 

We will consider in this article a family of stochastic processes which contains the 
classical squared Bessel processes. We recall that for 5 > 0, a 5-dimensional squared Bessel 
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process started at x > is defined as the unique strong solution of the stochastic differential 
equation 

dX u = 5du + 2 \J X u dW u and Xq = x > 

for all u € M + where (W u ) u >o denotes a Brownian motion. There has been a huge number 
of publications concerning this class of processes, resulting in many tractable tools for its 
study. We refer in particular to Shiga- Watanabe [14], Pitman- Yor [11], Revuz-Yor [13] and 
Deelstra-Debaen [4] , [5] . A natural extension of this family is to replace the dimension 5 by 
a function of the time variable. This was done by Carmona in [2], extending some of the 
previously mentioned tools to the class of squared Bessel processes with dimension 5 which 
is a time varying function. A second generalization of this class of processes was done by 
Deelstra and Delbaen in [3], [4] and [5]. They introduced a constant drift parameter (3 and 
used a dimension 5 which may be stochastic. Indeed, they extend results of existence of 
the solution of an extended stochastic differential equation 

dX u = (5 U + 2/3X u )du + 2y/xZdW u and X = x > 

They proved that the properties of scaling and additivity of the process (X u ) found by 
Pitman and Yor in [11] for classical squared Bessel processes stay valid for these extended 
squared Bessel processes. Moreover, Pitman and Yor in [11], using some classical space- 
time transformations and Girsanov's theorem, have shown that the study of squared Bessel 
processes with drift can be reduced to the study of the squared Bessel process without 
drift (i.e. j3 = 0). This property is equally true in the case of non-constant dimension. 
Finally, Delbaen and Shirakawa in [6] applied the class of squared Bessel processes with 
time varying dimension 5 to the square root interest rate model. 

Hence our contributions are: firstly, we generalize the family of squared Bessel processes 
to the class of squared Bessel processes with stochastic dimension (S u ) and with a drift 
parameter j3 u which is a function of the time. We will call this kind of process a generalized 
squared Bessel process (GBESQ) with stochastic dimension (S u ) and with a drift function 
f3 u , defined as the unique solution of the stochastic differential equation 

dX u = {5 U + 2/3 u X u )du + 2y/xZdW u and X = x > 0. (0.1) 

We establish several classical properties for this generalized process, extending or sometimes 
overlapping results from the aforementioned articles. This includes results of existence and 
uniqueness of the solution of (0.1); scaling and additivity properties of this solution. We 
finish by showing that, as in the two other previously mentioned settings, one can reduce 
to the situation (3 = through a change of law. Therefore we make a more precise study 
of the un-drifted case and give some rather tractable tools in order to compute the law of 
additive functionals of this process. 

We then give a more precise and theoretic version of the last result, in the form of a 
Levy-Ito representation of this process, similar to the one exposed in Pitman and Yor [11]. 
However due to the fact that the dimension parameter belongs to a function space and not, 
as in the classical setting, to M + , we cannot extend to our case all the results of global 
monotone coupling that exist in the classical setting. Our tools rely mostly, apart from the 
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probabilistic tool used in the setting of Bessel Processes, on functional analysis in order 
to extend results on the half line to a functional space. The Frechet-Riesz representation 
Theorem is a quite typical example of the theoretical results we use. 

Then we introduce the notion of Generalized Squared Bessel Bridges by conditioning 
our process to its endpoints. We give some extensions of the classical properties of Bessel 
bridges to our case, including an extension of the Bessel Bridge decomposition of Pitman 
and Yor [11]. We use this decomposition in order to give a general formula for the Laplace 
transform of additive functionals of the Bessel bridge. We give an example of how we can 
use this formula to get an explicit result. The techniques we use are very tractable and 
can be applied whenever one knows how to handle classical Bessel processes. We give a (if 
not very tractable) complete formula covering most settings in Appendix A. Our formula 
contains a power series, the value of which we are unable to compute explicitly, we however 
give some estimates on the speed of convergence of this series. 

We conclude this paper by giving some applications to financial mathematics of our 
results. Indeed, we will first give some examples of financial models which are GBESQ. 
Then, we will give three examples of applications: first we will give an explicit way to 
simulate stochastic volatility model where the volatility process is again a GBESQ, then 
the evaluation of zero coupon bond where the interest rate is a GBESQ, and finally, we will 
give an explicit method to simulate default times in credit risk model using a stochastic 
default intensity which is a GBESQ process. 

1 Squared Bessel processes with time varying dimension. 

1.1 Squared Bessel process 

Let (p,, (.7 : i)t£[o,T])ff > ) be a filtered probability space and let (Wt) t6 r 0) Ti be a Brownian 
motion under P. We suppose that the filtration {J T t)t&[Q,T\ ls generated by {W s ;0 < s < t}. 
We first recall some results of [13], about squared Bessel processes with constant dimension 
5 > 0. 

Definition 1.1. The squared Bessel process of dimension 5 > 0, denoted by BESQ% is 
defined as the unique strong solution of the stochastic differential equation (SDE) 

dX u = 5du + 2^fx~ u dW u and X = x > 0. (1.2) 

In this case X u > 0, for all u G [0, T]. 

We will denote in the sequel by Q s x , the law of BESQ^. We now recall the property of 
additivity of squared Bessel processes which can be found in Theorem 1.2 Chap XI of [13] 

Theorem 1.1. For every 5,5' > and x,x > 0, we have 

Q S X ®Q S '> = Q S+ J,. (1.3) 

X x+x v ' 
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1.2 Generalized squared Bessel process (GBESQ) 

Let C be the space of continuous functions M + — > M + and D the space of continuous 
by part, measurable and locally bounded functions M + — > R + . 

Definition 1.2. For all u E M + , 5, (3 G D, we will call the solution in C of the stochastic 
differential equation 

dX u = {2p u X u + S u ) du + 2^fx~ u dW u (1.4) 

with initial value X$ = x > 0, a generalized squared Bessel process with dimension (5 U ) and 
drift function (/3 U ). We will denote X as a GBESQx' 6 ■ 

We will denote by PQ S X the law of X and Q 5 X :=° Q s x . 

Remark 1.1. In the above definition (8 U ) and (/3 M ) are deterministic functions of the time. 
It is obviously possible to extend these definitions to random (5 U ) and ((3 U ) provided their 
trajectories belong almost surely to D, by conditioning on the realization of (5 U ) and (/3 U ). 
Thus all the following result are true conditionally on (5 U ) and (/3 U ), or, using the typology 
of random media, in the "quenched" setting. It is also straightforward by time restriction 
of the above SDE that the solution in this case is progressively measurable with respect to 
o-((d~ u ) u <t, {Pu)u<U (W u ) u <t). 

1.2.1 Existence, unicity and positivity of the solution 

We will obtain the uniqueness of the existence of the GBESQ and its positivity for a 
more general volatility structure. Indeed, we take the process X given by 

dX u = (2p u X u + 5 U ) du + a(X u )dW u (1.5) 

where a : ]R + — > M + is a function, vanishing at zero and satisfies the Holder condition, 
namely that for all x,y £ M + , there exists a constant c such that 

\a(x) - a(y)\ < cyj\x - y\. (1.6) 

Remark 1.2. 1. We can remark that the volatility structure of our GBESQ process 
(1.4) satisfies this condition. 

2. This class of processes was studied by Deelstra in [3] with a constant drift ft and a 
dimension 5 such that 5 u du < oo for all t £ M + . 

3. If 5 is in D then 8 satisfies point 2. 

Proposition 1.1. Let f3 and 5 be in D and a a function satisfying (1.6), then for every 
x > there is a unique solution Xt of the stochastic differential equation (1.5) such that 
X = x. 

Proof. The proof follows from Theorem 3.2.1 in [3]. Whereas as we said in Remark 1.2 
point 2. our model has a non constant drift term. But since f3 is in D, then it is locally 
bounded. Hence we can work with this new drift in the proof of Theorem 3.2.1 and obtain 
the same result. □ 
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To obtain the positivity of this solution, we will use a comparison result. 

Proposition 1.2. Let X 1 and X 2 be two solution of (1.5) with dimension 5 1 G D and 
5 2 G D and respectively initial state x 1 > and x 2 > 0. VFe suppose too, that the processes 
X 1 and X 2 are associated with the same Brownian motion (W u ) and the same drift (3 G D. 
If for all u > we have x 2 > x 1 and 5 2 > 8\ a.s., then for every u > 0, we have 

P [X 2 > Xl] = 1. 

Proof. We will again follow the proof of Theorem 3.2.2 in [3] with a non constant drift f3 
using the fact that f3 u is in D for all u G [0, T]. □ 

Corollary 1.1. The process X is always positive. 

Proof. This follows directly from the Proposition 1.2. Indeed, taking Xq = and 5 1 = 0, 
then the unique solution of (1.5) is given by X 1 = 0. Then taking X 2 another solution of 
(1.5) such that X$ > X^ = and <5£ > ^ = a.s. we obtain that X 2 > a.s. □ 

Remark 1.3. Hence taking for all u G [0,T], o~(X u ) = 2y/X u , then we have that the 
stochastic differential equation (1.4) admits a unique positive solution. 

1.3 Results about GBESQ 

We come back now to our GBESQ process given as the unique solution of the stochastic 
differential equation (1.4) 

Proposition 1.3. Let (Xu^) u ^[o,T] an ^ (^u ^)u€\o,T] be two independent generalized squared 
Bessel processes with time varying dimension S U ,5 U > with the same drift (3 U which is a 
bounded function and initial states x,x > 0. Then 

{JfW +y/);0 < u,X«> = x,Y^ = x} ^ < u,H^ =x + x}, (1.7) 

where H is another GBESQ with stochastic dimension 5 U + 8 U) with drift equal to (3 U and 
initial state x + x . In other words, 

Proof. We proceed similarly to the proof of Theorem 1.2 of [13] which is done for simple 
squared Bessel processes. Let W and W be two independent linear Brownian motion. Let 
X and Y defined by (1.4) as the corresponding two solutions for (x, (3, 5) and (x , (3, 5 ) and 
set H=X+Y. Then for any u G [0,T], 

f*U, f*\L (*\L 

H u = x + x + J 2/3, (X s + Y s ) ds + J (s a + 5'^jds + 2j (yx~ s dW s + y/Y,dwty . 
Let W be a third Brownian motion independent of W and W . The process C defined by 

r vx~ s dw s + VT s dw' s r „ 

Cu- 1{H S >0} jw 1 +/ l{H s =o}dW s 

JO V"s JO 
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is a linear Brownian motion since (CjC)« = u an d we have 

PU PU f'U 

H u = (x + x') + J 2p s {X s + Y s )ds + J + + 2 J ^/W s dCs 

f'U PU 

H u = (x + x) + J (2p s H s + (S s + <£) ) da + 2 J y/H~ a d£ a , 

which completes the proof. □ 

Proposition 1.4. For all u € [0,T], If X is a GBESQx' 6 (i.e. X = x), then for any 
c > 0, the process \X CU is again a GBESQ C £ ,S with 8 U = 5 CU and U = (3 CU . In terms of 

c 

distribution, this reformulates to 

Proof. We know that X u = x + J Q U 5 s ds + J™ 2j3 s X s ds + 2 J Q M yfX~ s dW s . So by a change of 
variable in this stochastic integral, we obtain 

u u u 

X u = x + c5 cv dv + / 2c(3 cv X cv dv + 2 \J X cv dW cv . 
Jo Jo Jo 



Hence 





X 






1 x - 

— si-u — 




+ 


/ 5 cv dv + 


C 


c 




to 




x 




u 


— si-u — 
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/ 5 cv dv + 


C 


c 




to 




X 




c u 






+ 


/ 5 cv dv + 


C 


c 




to 



-X cv —=i 

c vc 



-X cv ^-=dW cv . 

c vc 



Moreover, we know that if for all u £ [0,T], W u is a Brownian motion then, for all c > 0, 
W u := -±=W CU is again a Brownian motion. Hence this stochastic differential equation 



becomes 



^dX cu = (2(3 CU X CU + 5 CU ) du + 2^X cu dW u = \ 2cf3 c J-X cu + 5 cu j du + 2\^-X cu dW u . 

Consequently, the result follows from the uniqueness of the solution to this stochastic 
differential equation. □ 

1.3.1 Removal of the drift term. 

We now give a generalization of the change of law property which can be found in [11]. 
This proposition allows us to vanish the deterministic drift function. Let X be a GBESQx' S 
and for t G K + define the process Y by Y% = (3 s \/X s dW s , where the function /3 is in D 
and assume moreover that j5 is differentiable. Then the continuous local martingale defined 
by 

l 



Z t =£(Y t )=exp{Y t --(Y,Y) t 
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is equal to 



> J p s ^/T s dW s - p 2 s x s ds 



exp 



t 2p s ^Y s dW s + [ Msds- I Msds- [ (3 2 s X s ds 
o J o Jo Jo 

t rt rt 



P s 5 s ds + 2y/X 8 dW s 



p s 5 s ds- / pjX s ds 



Since X is a GBESQx' S , we obtain that 

1 



exp 



t rt 

(3 s dX s - / f3 s 5 s ds 
o Jo 



(1.9) 



By integrating by parts, we have 



ptX t = p X + f p s dX s + I X s dp s = /3 X + f p s dX s + I X s p' s ds. 
Jo Jo Jo Jo 

Hence by substituting this expression in (1.9) we obtain 



Z, 



expi- 



exp 



P t X t - /3 X - / X s (3 s ds - / p s 5 s ds - / piX s ds 
Jo Jo Jo 

PtX t - /3 X - J Msds - J (p' s + /3 S 2 ) X s ds 



If we show that Z± is a true-martingale then we will can apply the change of law formula 
given by the following proposition 



Proposition 1.5. 



Z t := 



Q8 " ex P 



t ft 

PtX t - /3 X - I Msds - ( a + l s ) X s cte 
o Jo 



(1.10) 



In fact, since X® ' S is positive for s G [0, T] and (3 G L>, so /3 is locally bounded. We have 
that ( /3 S + /3 2 ) > then the local-martingale Zt is bounded, therefore it is true-martingale. 



1.4 Some law results 

Starting from now and having in mind the previous result, we will always assume that 
the drift term is zero. We first recall a result from [2], Proposition 3.4 

Proposition 1.6. If 5 G D and n is a positive Radon measure on [0,a[, then 

Q s x exp (- £ X v dp(u)J = exp (| <&>)) exp (± £ §^<W«J , 
where is the (well and uniquely defined) solution in the distribution sense of 

= $n. (1.11) 

which is positive, non-increasing on [0,oo[ and follows the initial condition $^(0) = 1. 
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Remark 1.4. 1. With straightforward calculation, we can obtain a extended version of 
this result with non zero drift term, namely 



Pq 5 



exp 



X u d(j,(u) 



exp (f $») exp 11 J 



1 f a ^^n) 

d u au 



^ 



where the unique solution of the differential equation, in the distribution sense, 



d/J, u = i^-du + Pu-^du 



(1.12) 



which is positive, non-increasing on [0, oo[ and follows the initial condition $^(0) = 1. 



2. All these results could be extended under proper assumptions to measures without 
compact support. As our goal is to study the Bessel bridge, we feel free to do this 
restriction. 

The previous form gives a precise explicit expression for the Laplace transform of ad- 
ditive functional of the generalized squared Bessel. However in practical cases this explicit 
form is rather difficult to use. We give a more abstract form which, surprisingly, suits better 
our goals, and allows us to introduce some notations. 

Proposition 1.7. Let X s be a generalized squared Bessel process with dimension 5 and 
started in x. If 5 £ D and fi is a positive Radon measure with compact support on M + , then 
there exist a constant A^ and a function B^ : M + — > ]R + such that 



exp 



X u d/j,(u] 



A* exp 



Bfj,( u )d~ u du 



Proof. It is a trivial corollary of Proposition 1.6. 



□ 



Remark 1.5. Note that, by taking 6 constant, the A^ term is the same as in Corollary 1.3 
of page U0 of [13]. 

We can apply this formula to measures of the form \x t '■= Ae^, where et stands for the 
Dirac measure in t. Our strategy to compute A^ t and B^ t will be to take well-chosen test 
values and functions for x and 5. Taking 5 = 1, A^ t was computed in [13], page 441, and 
is equal to 



.4 



A 



(1.13) 



To get access to B^ t , we take x = and use test functions of the form 5 = X[b,oo] where 
X denotes the indicator function. For such 5 and x = 0, the generalized squared Bessel 
process is easily described: indeed it is identically equal to for < t < b and then evolutes 
like a squared Bessel process with dimension 1 started at 0, that is, the square of a standard 
Brownian motion. Therefore straightforward computations imply 



exp 



X u dfM 



En 



e~ XB l-»)+ 



1 



yi + 2A(£-6)+' 



(1.14) 
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(for the last equality see for example p441 of [13].) On the other hand this is also equal to 
exp B fj,(u)\. Taking the log and then the derivative with respect to b, then bringing 
everything together, we get the following result : 

Corollary 1.2. Under the statement of the Proposition 1.7, we have 

{-i^M-jfi+fc)*}- 

This is the Laplace transform of the transition Kernel of Q x . In the case of constant 5 
an inverse to this transform can be computed and therefore an expression of the transition 
Kernel, it is however not generally possible for any 5. We give in Remark 3.8 a simplification 
of this formula in the case of piecewise constant 6, and in section 5.1 some hints on how to 
overcome this issue in numerical simulations. 



2 The Levy-Ito Representation of GBESQ 

We now turn to the so-called Levy-Ito representation. For a motivation of such repre- 
sentation we refer to §4 of [11]. 

Theorem 2.2. 1. There exists a unique measure M on C such that for every random 
variable I on C of the form I = J Xtn(dt), with fi a positive Radon measure with 
compact support on (0, oo) and every a > 0, denoting by M u the distribution of 
(Xr t _ v \+ ,t > 0) for (Xt,t > 0) distributed according to M we get 




(2.15) 



Moreover M corresponds to the °M of Theorem 4-1 of [11]. 
2. For any 5 € D, there exists a R + — > C process (Y%,x > 0) such that 

Y X S =Y° + Y S =:Y X + Y\ (2.16) 

where (Y x ,x > 0) and (Y s ,5 £ D) are independent, Y x is as in Theorem 4-1- of [11] 
and Y s is distributed according to Qq. 

Remark 2.6. 1. This theorem should be understood as a Levy-Khinchine formula on 
C. For more on this we refer to §4 of [11]. We also refer to §3 of [11] for background 
on the measure M . In short, M is the excursion law associated to the zero diffusion 
Q x renormalized so that its entrance law is given by 

M(X t £dx) = ^e-f t . 

In particular the associated transition kernel starting for any strictly positive time is 
the same as the one of the squared Bessel process of dimension 0. 



9 



2. We also recall the expression ofY x from [11]. Namely, 

Y x (0) = x, Y x = J2^v on (0, oo) 



v<x 



where A v is a Poisson point process on C with Poisson measure M . Note that, due 
to the form of M , the sum in the last expression is actually infinite, however, for 
every finite t, the set of A„ such that A v (t) > is finite and its cardinality follows a 
Poisson distribution of parameter Jr. 

3. In [11], the representation in terms of a Poisson point process is also valid for S. 
However we do not have such a representation in our setting and believe it is a difficult 
question. It is however not needed for our purpose. 

Proof of Theorem 2.2: First, we would like to prove (2.15). It is clear that both terms are 
multiplicative in 5 and positive. The form in Proposition 1.6 also easily imply that the 
left hand side is continuous in 5 with respect to the infinite norm. The continuity of the 
right hand side follows easily once one knows that 5 u M u du (e~ aI — l) is finite, which 
is easily deduced from Theorem 4.1. of [11]. Hence we just need to prove it for well chosen 
test function. We take for test functions the functions that are constant on an interval and 
equal to zero outside. It is easy to see that such function form a generating family of D. 
Indeed it is rather easy to see that any constant by parts can be constructed as a sum of 
such functions. The fact that constant by part functions are dense in D is classical. 

More precisely, let 5 := dx[ a ,b]- We can easily see that for such 5, X t is identically 
zero on [0,a] then evolves like a d-dimensional squared Bessel process on [a, b], then like 
a 0-dimensional squared Bessel process. It is the clear that by a shifting of time we can 
reduce to the case a = 0. The Markov property implies: 

Q x [ e ~ aI ] = Q d x expj-ajf X t dfi(t)\ E° Xb expj-a^ X t dfi(t + b) 

We write I\ := J Xtdfi(t) and I2 = Jq°° Xtdfi(t + b). Then the previous expression can be 
written 

Qt [e~ aI ] = Qt [e- ah Q°x b [e~ ah ]] ■ (2-17) 
Theorem 4.1 of [11] implies that the previous expression is equal to 

gd L -a(I 1 +X b M(e- aI 2-l)) 

Denoting by /[ the function on C defined by 

I[((X t ) t > ) = Ji((% ) + X b M(e~ ah - 1) 
and applying once again Theorem 4.1. of [11], we get 

Q s x [e~ aI ] = exp I (xM + d M u du^ {e~ a ^ - 1) j . 
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It is clear from the previous form of I[ that for every u > b, M u (e aI>1 — 1) = 0, therefore 
the last expression can be reduced to 

Q s x [ e ~ aI ] = exp | (xM + d M u du^j (e _aJ i - 1)1 . 

We now want to prove that the above expression equals 

exp | (xM + dj^ M u di?J (e~ aI - 1)| . 

For this we only need to prove that for every < u < b, 

M u (e- aI 'i -l) = M u (e- aI - 1). 

We get this basically by doing the previous computations backward. Indeed 

r-6 



M u (e 



-ali 



1)) 



= M u 


exp < 


r\ 


/ 








/o 


= M u 


exp < 


r\ 


/ 








/o 


= M u 


exp < 




/ 











where we used in the second line Theorem 4.1 of [11] and in the third line the Markov 
property for M u and the fact that, by definition of M, for u < b, the transition kernel of 
M u for times greater than d is equal to the transition kernel of Q°. 

Now we need to show (2.16): it is trivial, just by taking Y x as in [11] and any independent 
5— squared Bessel process started at zero for Y s . Then Proposition 1.3 gives the result. □ 

3 Bessel bridge Decomposition 

We now turn to our main interest, namely the squared Bessel bridge with non-constant 
dimension. We define for 8 € D the law 



Q^Uy ■= Q d (-\X = x,X t = y). 



(3.18) 



In all the result we are going to mention, we are only interested in the trajectory up to time 
t G [0, T] of the Bessel bridge. So with a slight abuse of notation we will consider Qx^-y as 
a probability measure on C([0, t], R + ). We will also make use of the time reversed Bessel 
bridge, namely we call Q the Q distribution of (X(t — s))o< s <t- In particular it is easily 
seen that we have 

Lemma 3.1. Let X be a GBESQ S,X such that for t > 0, X t = y, then we have that 

A5,t g5,t 



where 5{s) = 5(t — s). 
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We also restrict to the case where 5 is not equal to zero excepted on a discrete set of 
points, as it will be the case in financial applications. This is only to ensure that for every 
y > 0, the conditioning event in (3.18) has zero probability. It will also ensure that the 
distribution Q x 'Uy is continuous with respect to (x,y). 

We recall that, by definition, Qq'X v is the time reversal of Qy'^o- We first give a "bridge" 
version of the additivity property. 



Proposition 3.8. Let 5, 5' be in D, then 

Q S,t Q 5',t _ Q S+S',t 

Remark 3.7. In this setting there is no straightforward representation in term of the 
unconditioned process in the sense of (5. a) of [11], whereas doubtless a careful application 
of the arguments of [18] should lead to some similar (but not equal) result. This is however 
not needed for the above proposition to be true. 

Proof of Proposition 3.8. This is directly obtained by conditioning the statement of Propo- 
sition 1.3, using the (obvious) fact that, for X and Y two M + -valued processes, the events 
{Y(t) + X(t) = 0} and {X(t) = 0, Y(t) = 0} are equal. □ 

We are aiming at obtaining an extended Bessel bridge decomposition. Firstly, we need 
some notation and an introductory Lemma. We denote by the inverse Laplace trans- 

form of /. And we call 

Fg(\) = exp ( ^—) exp (- f ^ -du) . (3.19) 



1 + 2Xt J r { J 1 + 2A(t - u) 



Lemma 3.2. Let 

^-£- X n + 2n(y) 



(x/(2t)r r -l Ft 

b^Jn 



5 > X > y{ ' ' W(2i)) fc /.-] F t t,.\ 

then bfixyfo) is the distribution of U given 

X*+X 1 + --- + X u = y, 

where U ~ Poisson(x/ (2t)) , X\ ~ exp(l/(2t)) are i.i.d. random variables and X* is a 
random variable independent ofU and Xi and its Laplace transform is given by F$(\). 

Proof. By simple computation one gets that the Laplace transform of X\ is equal to 

expi— f — -du 



o l + 2A(t-u) J 1 + 2U 

We deduce that, conditionally on U = n, the Laplace transform of X* + X\ + • • • + Xjj is 
equal to Fj +2n {X)- The continuous version of the Bayes Formula then gives the result. □ 

Then we obtain the extended Bessel bridge decomposition 
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Theorem 3.3. For every 5eD,xeR + ,yeR + , 

oo 

Qx— >y 



n=0 



where the last term on the right is the mixture of the laws Qq^q, n = {0, 1, . . . } with weights 
given by b\ x y (n). 

Proof. Using point 2. of Remark 2.6, we get that, for Y* ~ Q s x , 



Y X 5 = Y X ° + Y S , 



where the two processes are independent, Yq ~ Qq and 

Y° - VA 



v<x 

where A v is a Poisson point process with measure M. We separate as in [11] this sum in 
two, keeping on one side the terms which are zero at t and on the other side the terms 
which are positive at t, we get that 

x x — x x0 "r i x + ' I i 

where the three processes are independent, 

Y x0 ~ ^^0' 

and 

where U is a Poisson(x/(2t)) variable independent of everything else and Zi are i.i.d. 
trajectories whose distribution is given by M conditioned on Xt > and renormalized. 
We condition on the values of Zi and Yq at time t. Conditionally on U = u, and Zi(t) = 
Zi, Vi < U, it is shown in [11] that 

Summing this, we get that, conditionally on U = u, and Zi(t) = Zi, Vi <U with ^ Z{ = z, 

Y x+ ~ Q _> z Q ^ . 

On the other hand conditionally on U = u, and Zi(t) = Zi, \/i <U, using time reversal and 
Proposition 3.8, 

Y s ~ Q 5 / = Q 5 / n ffi Pn'* 
Putting everything together, we get that conditionally on U, 

n&,t _ Q$,t ~ Q0,t m Q0,t ffi n 4 "'* 

On the other hand, conditionally on Y, Y^t) has Laplace transform FUX) and Zi(t) are 
independent random variable with law exp{l/(2t)} Therefore, using Lemma 3.2, we get 
the result. 

□ 
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Remark 3.8. The expression in Lemma 3.2 is fairly general, however in most practical 
applications 5 is constant by parts. We thus give a simpler expression in this case. Indeed, 
let = ro < T\ ■ ■ ■ < r n = t be an increasing sequence of times such that on each [7$, Tj.fi), 
5 is constant equal to Si, then 



*i=->{-T£»}ll 



Xx Ifl ■( l + 2A(t-T i+ i) 
1 + 2A(t - Ti] 



3.1 About the coefficients in the Bessel bridge decomposition 

The main difference between the classical setting and ours is the fact that the coefficients 
^8xy( n ) are n °t as explicit. However they can be computed for given values of 5. For 
numerical purposes it is interesting to have a bound on the speed of decay of these terms. 
We have the following result: 

Proposition 3.9. Let 5 G D, d = inf^t] 5 andx,y £ ]R + , then for every increasing function 
f : N h-> M + , and with u = d/2 — 1 and z = J~xy, 



E *W*)/(*0 ^ E ^Ak)f(k), (3.20) 

where 



MQ = (|) 



r(n + z/ + 1)^(2)71! 

are i/te coefficients of the Bessel(v, z) distribution (I v {z) stands for the modified Bessel 
function). In particular the power series associated to the bg x (k) has an infinite radius of 
convergence. 

Remark 3.9. The philosophy of this proposition is that, when one wants to compute quan- 
tities related to the Bessel bridge, the truncation arguments used in the classical setting are 
also valid in our setting. 

Proof. From the representation of Lemma 3.2, it is clear that bg x (k) is the distribution 
of a variable U which is negatively correlated with X*, the distribution of the generalized 
squared Bessel process at time t. Proposition 1.3 implies that we can couple the generalized 
squared Bessel process of dimension 5 with a classical squared Bessel process of dimension 
d, the first one being always bigger than the second one. From this we deduce that the 
distribution of U for 5 is dominated by the distribution for d. This distribution is computed 
in [11], (5.j) and is the Bessel(u, z) distribution. The result follows easily. In particular 
the parameter t does not intervene, as it simplifies between a and A in [11], (5.j). □ 



4 Laplace transform formula associated to extended Bessel 
bridges. 

We are now interested in computing the Laplace transform of some functionals of the 
Brownian bridge. We first deal with the simpler case of the bridge to zero. In this case, 
because additivity holds, things are simpler. We have the following: 
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Proposition 4.10. For every Radon measure [i on [0,t], there exist Aq{o) and a function 
Bq(u) : [0,t] —> R depending on fi such that 



Q 



8 A 



exp 



aX u dfi(u 



(A (a)) x exp{ / B%{u)5 u du 



Remark 4.10. Note that, for t = 1, Aq{ol) is as in in (5.d') of [11]. We will give an 
example of explicit computation of Aq(o) and Bq{u) in Theorem 4-4- 



Proof. We first assume that the function $ defined by (1.11) is such that j Q 
Using Proposition 3.8, it is easy to check that 



t ( 



du < oo. 



(x,5) :-- 



log Q^o 



exp 



aX u d/j,(u) 



is additive in both terms, therefore (f)(x,5) = 4>(x,0) + 4>(0,5), and 
is necessarily of the form 4>(x, 0) = oqx. We now want to show that 



(x, 0) being linear, it 



log Qq'^o 



exp 



aX u dfi(u] 



is of the form "J™ BQ(u)5 u du" . To show this we need to apply the Frechet-Riesz represen- 
tation, which states (in our context) that any continuous linear form on the Hilbert space 
L 2 ([0,i]) is of the above mentioned form. So, keeping in mind that D is a dense subset of 
L 2 ([0,i]), the only thing we need is to show that <p(0,5) is continuous with respect to the 
L 2 norm on [0, £]. As we now it is linear, it suffices to show that it is bounded on the unit 
ball. By proposition 1.6, and the Cauchy- Schwartz inequality, we now that it is true for 
the unconditioned process, as by assumption ^ is square integrable. Therefore the result 
holds as a consequence of the following coupling result, whose interest goes further than 
the context of this proof. 



Lemma 4.3. For every 5 G D, x S 
such that 



there exist a 



') valued process (X S ,X® 



s )s€[0,t] 



x x 



Q 5 X° ~ Q 

■<-X1 s -< 



J^ , and X Q S < X s a.s., V0 < s < t. 



Before giving the proof of this lemma we finish the proof of Proposition 4.10. We only 

rt /d>'\ 2 

1 du < oo. First note that by definition of <3?, for 



need to get rid of the assumption J| 



o V * 



all e, e (^-) du < oo. We deduce from the same arguments as above that for all e, there 



exist B Q ' such that 



a 



6j 



exp 



t-e 



aX u dfi(u 



(A (a)) x exp 



t-e 



Bq' € (u)5 u du 



From this it is clear that for e' < e, the restriction to [e,t — e] of B ' e is -B ' e - We- can 
accordingly define Bq on (0, t) such that the restriction to [e,t — e] of Bq is Bq' 6 . One can 
then easily check that Bq verifies the statement of Proposition 4.10. □ 



We now turn to the Proof of Lemma 4.3. 
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Proof of Lemma 4-3. The proof follows a standard scheme, see for example pages 205-206 
of [7] . The idea is to express the conditioned process as a solution of a SDE whose terms are 
smaller than the SDE defining the unconditioned process, and then apply the Comparison 
Theorem of [17] to prove the existence of a coupling. 

We have first that for s < t, denoting by qt(x, y) the transition kernel associated to Q s x . 
Using 1.3, we can write 

q t (x,y)= (q?(x,.)*q S t (0,.))(y), (4.21) 



where * stands for the convolution, and q® is the transition kernel of the 0-dimensional 
squared Bessel Process. In particular 



We have 



dQ 



s,t 

x->0 



dQl 



q t (x,0) = q°(x,0)q 5 t (0,0). 
<h-.(X.,0) 



(4.22) 



Qt(x,0) 



exp (logq t - s (X s ,0) - \ogq t (x,0)) . 



We write h(x,s) = \ogqt- s (x, 0). We now from (4.22) that 

h(x,s) = h°(x,s) + h 5 (0,s), 



(4.23) 



where h s (x,s) corresponds to the squared Bessel process of dimension 5. 

From this form it is obvious that the derivative of h with respect to x does not depend 
on 5, therefore is equal to its value for the 0— dimensional squared Bessel process, that is 
2(s_ t -) (see [7] for details), and in particular the second derivative with respect to x is . 
We have thus, applying Ito's formula that 



dQ 



s,t 



exp 



1 f s dh 

— rdX u + / —(X u ,u)du 

o 2{t - u) J Q ds 



Using once more (4.23) we see that 4^ can be decomposed into a part which only depends 



ds 

on 5 and another one which only depends on x. Putting everything together we get 

/ T 1 f s X" . f s dh s 

= eXP (-/ n 2ii 



dQi 



o 2{t-u 



■dX u 



, 2(t- u y du + J ^ u)du] 



We do not need to compute the last term in the exponential, indeed it is a deterministic 
term, and we now that the Radon derivative is a martingale, therefore its value is imposed 
by this condition. Recalling that by definition 

dX u = 5 u du + 2sfx~ u dW u , 

and replacing in the above expression, we have 



dQx^o 



dQl 



exp 



X u 



(t-u) 



dW u ,+ 



Xn, 



o 2(t 



;du + a deterministic function of s 



u 
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As X s — Jq 5 u du is a martingale (under 0,%), and keeping in mind that the deterministic 
term is imposed by the martingale condition, we have: 



dQ 



s,t 



X u 



dW,, 



o (t-u) 



Girsanov's Theorem then implies that, under Q^I+q, 



dX u = 5 u du + 2JX u dW u 



2X U 
t — u 



du. 



Therefore we obtain the result by an application of the Comparison Theorem of [17]. □ 

We deduce from Proposition 4.10 and Theorem 3.3 the following 

Corollary 4.3. Let for fi a Radon measure on [0, t], I = f Xdfi, Aq(o), Bq(u) : M + — > E 
as in Proposition 4-10, and Aq{q) be the Aq{q) for the image of (i under the map s —> (t—s), 
let also Bo (a) = exp | J" * Bq (u)<iu| then we have 

Q^U y (e- aI ) = A (aTA (a)y expjjf' B%(u)5 u du} ^ ^ >3/ (n)B (a) 4 ^ . 



We now give explicit computations when [i is the Lebesgues measure on [0, t\. 
Theorem 4.4. Let X be a GBESQ S X and let /3 = V2at, then we have 



Q S l+ y (e- a ti x ° ds ) = exp {^±^(1 - /3coth/3) 

coth(v / 2a(t - u))^ 5 u du^j &$,*,,,(«) ^ 



o \2(t-u) 2 



sinh/3 



2ir 



Proof Here [i is invariant under the map s — > (t — s), so that Aq{q) = Aq{o). Corollary 
3.3 on page 465 of [13] gives for t = 1, b € K + and constant 5 that 



Q 



5,1 



-y J X s ds 



b \ 2 fx 
— — ) exp ( — (1 — 6coth b) 
sinhfe/ F V2 V ; 



Taking first a := j, then although this expression is given for t = 1 we can extend it to 
[0, t] using the fact that, for 8 constant, 



Q^Je-«foXM=Q^( e -«<lo 1 x. 



i ds 



Q 5,l ( -od*fZx.d3 



where first, we made a change of variable and then used the scaling of Bessel processes. 
Hence taking y = and 5 = 0, we obtain the expression of Aq{o) 



Aq{o) = exp — jl — y/2at coth \/2crf j . 



(4.24) 
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And taking y = and x = we obtain the expression of Bo (a) 



B (a) = J (4-25) 
y sinnv2ac 

In order to compute Bq(u), we take as before test functions for 5, namely we assume 
6 = X[a,tb w ith < s < t. We have by the same arguments as before 

By scaling, this is equal to Bj(a(t — s) 2 ), where Bg(a) is the value of Bo(a) for t = 1. On 
the other hand it is equal to exp (J"* Bff (u)du). Then straightforward computations give 

B%{u) = - ^ coth(^(t - u)). 

This, together with the previous computations, gives the result. Note that Bq(u) is not 
defined for u = t, however it is clearly integrable on [0,t]. □ 

Remark 4.11. The previous computations may be easily applied to get formulas for general 
additive functionals of the GEBSQ. We skip for the moment these tedious computations, 
but give the complete result in Appendix A 



5 Some applications 

As we said in the introduction, our work was motivated by applications to Finance. 
Indeed, our class of GBESQ contains various financial models. Hence we can give an 
answer to many financial mathematical problems using our Bessel bridge decomposition 
and Laplace transform formula results applied to particular GBESQ. We begin thus by 
giving some financial models which are in the class of GBESQ. 

Ornstein-Uhlenbeck (OU): Let (Vt)te[o,T] following an Ornstein-Uhlenbeck model given 
by 

dV t = liVtdt + adWt with V = x. (5.26) 
Then we have the classical representation in term of GBESQ given by 



X i J?- (e- 2 ^ - r 



V t = ef A - 

and X is a GBESQ 5,X with dimension 5 = 1 and initial value x = x 2 . 



(5.27) 



Cox Ingersoll Ross (CIR): Let (Vt^g^r] be a CIR process then it is the solution of the 
stochastic differential equation 

dV t = (q - pv t )dt + a^VtdWt with V = x. 
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where a E M + , f3 £ K, a > 0. Then we have by a time space transformation the 
classical representation in term of GBESQ given by 

Vt = e -^ X (^(e^-l)) (5.28) 



V4/3 

and X is a GBESQ 5 > X such that 

dXt = <Sd< + 2^/X t dW t 
with a constant dimension S = ^§ and W" (e' 3 * — l)^ = J * ^e^dW s . 

Constant elasticity of variance model (CEV): Let (V^g^T] following a Constant elas- 
ticity of variance model (CEV) given by 

dV t = (iVtdt + aVfdWt with < p < 1 and V = x. (5.29) 

where fi represents the instantaneous mean, aVf the instantaneous variance of V and 
p is called the elasticity factor. In the limiting case, p = 1, the CEV model returns 
to the classical Black and Scholes model. In the case where p = 0, it is the Ornstein- 
Uhlenbeck model and if p = i, it is a Cox-Ingersoll-Ross (CIR) model without mean 
reverting. Then we obtain that 

r n n _ 1 w2 , x \ 1 27 1 



2(1-P) 

and V is a GBESQ S,X with dimension <5 = 2 pZ\ an d initial value x = x~ 2 ^ p ~ l \ 

Extended CIR: Let (Yt)t£[o,T} following an extended Cox Ingersoll Ross model which is 
given by the solution of the stochastic differential equation (SDE) 



dV t = {a t - PtV t )dt + atVVtdWt with V = x. (5.31) 

where at, Pt > and at > are time- varying functions and at is continuously differ- 
entiable with respect to t € [0, T]. Now, we recall Corollary 3.1 of [15]. 

Corollary 5.4. We have 

V t = exp (- p u dv)j X Q jT *^dv\ , (5.32) 

where X is a GBESQ S ' X with time varying dimension St = with v t := 

2 Jo T^du, 0t = exp ^— Jq /3 u duJ and initial value x := = x. 

Hence, the class of GBESQ contains a large class of standard (financial) models 

Remark 5.12. We see from these example that one needs in many occasions to evaluate 
the law of junctionals of the form f g(s)Xfr s \ds for various functions g and f (with f 
increasing), for X being a GEBSQ or a Generalized Bessel bridge. Laplace transform of 
these quantities can be computed using the results of Remark 1.4 or alternatively Appendix 
A, using the simple argument that, by change of variable, 
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5.1 Generating values of X t given Xo and its Laplace transform. 

Let X be a GBESQ 5 > X . Our aim is to obtain a sample of generating values of Xt 
for t G [0, T] assuming that with know Xq = x and that we have an explicit formula of 
the Laplace transform of Xt- A classic way to do this is to inverse the Laplace transform 
formula. But in our case the Laplace transform formula we obtained in Corollary 1.2 is not 
easy to inverse. So we need an other method to generate a sample of values of Xt- One can 
use for example the method developed by Ridout in [12] which allows to generate sample of 
random numbers from a distribution specified by its Laplace transform (without having to 
inverse the Laplace transform). Hence, as we can evaluate explicitly the Laplace transform 
of X assuming Xq = x with our Corollary 1.2 and simplification given by Remark 3.8, one 
can apply the method of Ridout in [12] to obtain a sample of values of Xt for any t G [0, T]. 

This method also permits to jointly simulate Xt (for a fixed t) and any additive func- 
tional of X s , < s < t without simulating the whole process. To do this first simulate Xt 
through its Laplace transform, then plug this value into the formula of Appendix A and 
repeat the process. We give an example in Section 5.2 

5.2 Simulation of integrate process conditioning to its starting and end- 
ing points. 

In this application, we follow results found by Broadie and Kaya in [1] and extended 
by Glasserman and Kim in [8]. We will apply the same idea but using our more general 
framework. Indeed, assume that we would like to simulate a stochastic volatility model 
given for all * G [0, T] , 6 6 D by 



dS t = mStdt + sJVt [pdWl + Vl - P 2 dWn , S = s, (5.33) 

dV t = S t dt + 2^/VtdW} , V = x. (5.34) 

where W 1 and W 2 are two independent standard Brownian motion, Xt is a Markov jump 
process on finite space S := {1,2, ... ,K} and p G [—1,1]. So the stochastic volatility 
process is a GBESQ S,X . Let for all t G [0,T] the solution of the first stochastic differential 
equation in (5.33): 

S t = S exp (J p s ds -lj*V s ds + p^ VV s dW} + v 7 ! - P 2 jT Vv s dW^j . 

Moreover J Q * \/V~ s dW]. = | (Vt — Vq — J ' <5 s <is^ . Then we have that log is conditionally 
normal given J * V s ds, Vt and Vq. Hence 



log [ £ ) ~ M (Y Psds - l -J^ V s ds - P -[V t 



V - J S s ds\ , (1 - p 2 ) 2 J V s ds\ 



Moreover we can sample Vt using again method stated in subsection 5.1. Hence we can 
obtain a sample of couple (Vo,Vt). Then the problem of exact simulation of our model 
(5.33) is reduced to the problem of sampling from 

Wsds^Vt 
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Since V is a GBESQ S X , we know by Theor em 3.3 that for every 5 G D, x G K + , y G R + , 
there exist processes Xi, X2, X3 and X4 such that 

jf* u s( i s | v = x, v t = y) = x 1 + x 2 + x 3 + x 4 , (5.35) 

where Zi ~ Q°'^ , A 2 ~ Q° 'X y , X 3 ~ Qj^ and A 4 ~ £^° =0 ^^(n)^. Moreover, 
depending on the GBESQ V chosen, we can use Theorem 4.4 or Remark 5.12 and Propo- 
sition 3.9 to have a way of sampling from ( J * V s d,s|Vb, Vtj ; and so an exact simulation of 
our stochastic model (5.33). 

5.3 Evaluation of Bond price in a credit notation model with regime 
switching 

Assume that we are under an equivalent risk neutral probability P, then from the 
general asset pricing theory, we have that the discount bond price B(T) at time with 
maturity T > of an interest rate r is given by 



B(T) = E 



exp - / r s ds \r 



(5.36) 



Let the process r be given by the solution of the stochastic differential equation given by 
(5.31), we assume that the value of the parameters (at, at) is function for all t G [0, T] of the 
value of an exogenous continuous time on finite state Markov chain (S)t & \o j<\. Moreover, 
we assume that = 0. Hence the Markov chain 5 takes value in S := {1, 2, . . . , N} and 
could represent the credit notation of a firm A given by a credit notation entity. Then 
the process r could be the spread of the firm A and we would like to evaluate the price at 
time t G [0, T] of a bond of this firm with maturity T > 0. (See Goutte and Ngoupeyou in 
[10] for more details about this defaultable regime switching modeling). Hence our model 
becomes 

dr t = a(S t )dt + a(S t )VndW t , r = r (5.37) 
and so a, f3 and a are piecewise continuous functions. 

Remark 5.13. 1. The fact that we assume that S is a exogenous Markov chain implies 
that S is independent of the Brownian motion W for all t G [0, T\ . 

2. Assuming that S is exogenous has a economic sense, indeed in representation (5.37) 
there is two source of randomness, firstly the market risk modeled by the Brownian 
motion W and secondly an external expertise risk factor S given by an external entity 
of the firm A. This could be the credit notation representing the wealth at time t G 
[0, T] of the firm A or an indicator of the wealth of the global international economy 
at time t. 

Let us denote by S the historical values of S, i.e. S := (St;t G [0, T]). Using (5.32), we 
have that 



1 /•* 



( \ I °lds 



o 
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where X is a GBESQ S,X with time varying dimension Sf = with v t := | J* a^ds 

and initial value x = r. Hence we have two different ways to evaluate 



B A (T) = E 



T 



exp — / r s ds \tq = r; S 



o 



Firstly, we can use the Proposition 1.6 and solve the corresponding equation (1.11). Or we 
can use our Bessel Bridge decomposition to jointly sample the final value of the process X 
and the Laplace transform. Indeed, by Remark 5.12, we have in this model that 

1 /"* 

f(t) = u t := - J a 2 s ds and r t = X f ^ (5.38) 

and g(.) = 1. So we can apply Theorem A. 6, which is the extended version of our Theo- 
rem 4.4, to evaluate the Laplace transform of J^Xf^ds using our extended Bessel bridge 
decomposition where we take for Radon measure fj, the measure given for all i > by: 

Mt) = jy^dt- 

Hence to apply this result, we need to have a sample of value of = Xf(T) with respect 
to the starting point tq = Xq := x and the law of the associated GBESQ for r. For this 
we will use the method stated in subsection 5.1 since we have a Laplace transform formula 
expression given in Corollary 1.2. Moreover the quantity b^ x An) given by Lemma 3.2 is in 
this case simpler since we have Remark 3.8. Finally these allow us to obtain an evaluation 
of the Bond price B A {T) at time zero of the firm A using Monte Carlo. 



5.4 Simulation of the default time in credit risk model 

We work in the model of construction of Cox processes with a given stochastic intensity 
A. Hence we denote by A a stochastic default intensity of a firm or country A. We assume for 
all t G [0, T] that J Q \ s ds < oo and that J °° X s ds = oo. Then the default time r of the firm 
or country A is a real random variable which is given by the first time when the increasing 
process X s ds is below a uniform random variable U independent with A. Hence we get 

r := inf jt > 0; exp (- J X s ds \ < u\ . (5.39) 

So to simulate the default time requires to simulate the stochastic default intensity A 
together with its time integral. Hence assuming that A is a GBESQ S ' X , we can use again 
our Bessel bridge decomposition given in Theorem 3.3, our Laplace transform formula given 
in Corollary 1.2 and the method stated in subsection 5.1, to do this. 



A Laplace transform of additive functionnals of the GEBSQ 

In this appendix we give an explicit formula for the Laplace transform of additive 
functionals of the Generalized Bessel bridge. We first recall Corollary 4.3: 
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Corollary 4.3. Let for \x a Radon measure on [0, t], 1^ = J Xdfi, Aq{o>), Bq(u) : M + — > K 
as in Proposition 4-10, and Aq{ol) be the Aq(o) for the image of (i under the map s — > (t—s), 
let also Bo (a) = exp | J Q * Bq (u)c?u| then we have 

S^(e-^) = M<*) X M<*) V expjjT* B%{u)5 u du} [f^ ^ >a; >)B (a) 4 ^ . 

The idea, used before, is to take test functions for 5 to get the explicit value of Ao(a) 
and Bq(u). We recall the result in the case of Standard Bessel Process, as given for 
example on page 465 of [13], in the proof of Theorem 3.2., taking y = and recalling that 

iinx^o 7^7 ~~ y ' 

Theorem A. 5. In the case of constant dimension 5, we have 

\ <V 2 



~ 2 ds 



where, as usual, <f)n is the unique solution of (J)" = n<p which is positive, non increasing and 
such that (f>n(0) = 1; and = 4-. 

We first extend this formula to bridges of arbitrary length (keeping for the moment S 
constant). Let [i be a positive Radon measure on [0, t]. We have by change of variable, 
then scaling, 

By Theorem A. 5, we deduce that 

where ^ is the unique solution of </>" = t 2 fi(t-)(j). It is clear that = <^(i-), where $ is 
the unique solution of <fi" = fi<p which is positive, non increasing and such that 4^(0) = 1. 
Expressing the previous result in term of (ft we have, with obvious notations, 

^[^1 - {i (*<» + \ - J^) *} (j^) ^ 

We deduce from this expression that Aq(1) = exp | \ M^(0) + j — — t ^_ 2rfjj ^ | • We also 
deduce that Ai(l) = exp < \ I -F'(O) + 7 — ri -. 4 „ — ) > , where <£* is the unique solution 

{ 2 \ * 1 Jo <^( s ) ds J J r 1 

* f (fit 

of 4>" = n(t — ■)(/) which is positive, non increasing and such that ^(0) = 1; and F^ = 

We now turn to the computation of Bq(u). For this take 5 = X[b,t] - We have as previously 
explained, 
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exp<j / B%{u)du \ = Ql'Xo 



where 4> t f ] b is the unique solution of <f>" = fi(b + -)4> which is positive, non increasing and 
such that (/>^ b (0) = 1. therefore = ^jqjjy ■ Therefore we have 

-p{jf )] = (_|Z1_) ,/2 . (A , 0) 

We take the logarithm and differentiate with respect to b to get: 

BMu) = — r+F,u) ^ 1 



■ J*" 6 X s d»(s+b) 



t-b 



1/2 



Jo b( fy b ( s ) 2ds , 



Finally, using (A. 40), we have 



B (l) = expj^ Bl(u)du 



1/2 



Finally, putting everything together, we have: 
Theorem A. 6. 

~ exp £ (,'(„, + I - J ^ e Xp | (..(„, + i - 

i . „,., > *!,(«)- 2 ^ , , 1 l^„. , J t V"' 



cxp 



[ I , 1 ; — — w4 



where (fr^ is the unique solution of <fi" = ficj) which is positive, non increasing and such that 
Yu(®) = 1' 0u * s ^ e unique solution of <fi" = /u(i — -)4> which is positive, non increasing and 

such that A(0) = 1; and F*. = % 
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